Chiral susceptibility and the scalar Ward identity 
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The chiral susceptibility is given by the scalar vacuum polarisation at zero total momentum. This 
follows directly from the expression for the vacuum quark condensate so long as a nonperturbative 
symmetry preserving truncation scheme is employed. For QCD in-vacuum the susceptibility can 
rigorously be defined via a Pauli-Villars regularisation procedure. Owing to the scalar Ward identity, 
irrespective of the form or Ansatz for the kernel of the gap equation, the consistent scalar vertex at 
zero total momentum can automatically be obtained and hence the consistent susceptibility. This 
enables calculation of the chiral susceptibility for markedly different vertex Ansdtze. For the two 
cases considered, the results were consistent and the minor quantitative differences easily understood. 
The susceptibility can be used to demarcate the domain of coupling strength within a theory upon 
which chiral symmetry is dynamically broken. Degenerate massless scalar and pseudoscalar bound- 
states appear at the critical coupling for dynamical chiral symmetry breaking. 
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I. INTRODUCTION 

The analysis of colour-singlet current-current correla- 
tors 



(«(aOr«(a:) 9 (0)rg(0)), 



(1) 



where T is a Dirac matrix, or, equivalently, of the as- 
sociated vacuum polarisations, plays an important role 
in QCD because these quantities are directly related to 
observables. The vector vacuum polarisation, e.g., cou- 
ples to real and virtual photons. It is thus basic to the 
analysis and understanding of e + e~ — > hadrons and a de- 
termination of the strong running coupling on the pertur- 
bative domain [l], Q • Nonperturbative information is also 
available. For example, the mass of a hadron can be esti- 
mated in lattice-regularised QCD by analysing the large 
Euclidean-time behaviour of a carefully chosen correlator 
[H, lij . Correlators are also amenable to analysis via the 
operator product expansion and are therefore fundamen- 
tal in the application of QCD sum rules [5j]. 

A given vacuum polarisation is a function of the to- 
tal momentum, P, which is conjugate to the spacetime 
separation between the two currents. Its value at P = 
yields a vacuum susceptibility; namely, a measure of the 
response of the theory's ground state to a fluctuation in 
some external parameter. Herein we consider the chiral 
susceptibility, which in the chiral limit measures the re- 
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sponse of the ground state to a fluctuation in the current- 
quark mass. In the theory of phase transitions this is 
analogous to considering the response of a magnetisation 
to an infinitesimal external magnetic field. It is funda- 
mental to understanding the phases that may be realised 
in a given theory. 

In QCD the quark-parton acquires a momentum- 
dependent mass function, which at infrared momenta is 
~ 100-times larger than the current-quark mass. The 
Dyson-Schwinger equations (DSEs) 0] explain that 
this effect owes primarily to a dense cloud of gluons 
that clothes a low-momentum quark [8J, |9j . This marked 
momentum-dependence of the dressed-quark mass func- 
tion is one manifestation of dynamical chiral symmetry 
breaking (DCSB), which is the single most important 
mass generating mechanism for light-quark hadrons; e.g., 
it is responsible for roughly 98% of a proton's mass. 

The behaviour of the chiral susceptibility can be used 
to determine under which conditions chiral symmetry is 
dynamically broken; e.g., to demarcate the domain of 
coupling strengths upon which the phenomenon occurs. 
It can also be employed to analyse the nonzero tem- 
perature behaviour of QCD, and to determine the criti- 
cal temperature for chiral symmetry restoration and the 
critical exponents associated with that transition; e.g., 

Refe.Ha.lSlSQ- 

Our aim herein is to explain the essence and uses of 
the chiral susceptibility in QCD. In this connection it is 
important to note that when working with a vacuum po- 
larisation or susceptibility the question of an appropriate 
regularisation scheme arises [1J, [l5|, [l6| • It is moot for 
QCD in-medium because no new divergences are encoun- 
tered in a theory that is properly regularised in-vacuum. 
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In Sec. [II] we therefore discuss the question of regulari- 
sation for the in-vacuum theory and present a number 
of model-independent results. The scalar Ward identity 
plays an important role. In order to illustrate the appli- 
cation of the chiral susceptibility we employ two simple 
models for the gap equation's kernel. They are described 
in Sec. IIHI The difference is expressed in the form of the 
quark-gluon vertex and in Sec. lIVI we show that this dif- 
ference has only a minor quantitative effect. That both 
Ansatze can be employed self-consistently owes again to 
the scalar Ward identity. We wrap-up in Sec.lVl 



II. CHIRAL SUSCEPTIBILITY 

In QCD the vacuum quark condensate associated with 
a quark of flavour / can be written [l7| 



a / (m / ;C,A) = Z 4 (C,A)iV c tr D / S f (q;m f ;0, (2) 

Jq 

where: mf(Q) is the renormalised current-quark mass, 
with C the renormalisation scale; Z 4 is the Lagrangian 
mass-term renormalisation constant, which depends im- 
plicitly on the gauge parameter; and :— J A cPqj (27r) 4 
represents a suitable regularisation of the integral, with 
A the regularisation mass-scale, which is taken to infinity 
as the last step in a complete calculation. 

The dressed /-quark propagator, Sf in Eq. ((2|), is ob- 
tained from the gap equation 1 



Sfip)- 1 =Z 2 (i 1 -p + mf 1 ) + E/(p) , (3) 



with 



E/(P) = Zi f q 3 2D ^(P ~ ?)y7^/(«) yrSfap), (4) 

where D fll/ (k) is the dressed-gluon propagator, T^(q,p) 
is the dressed-quark-gluon vertex, and m^ m is the A- 
dependent current-quark bare mass. The quark-gluon- 
vertex and quark wave-function renormalisation con- 
stants, Zi^(( 2 , A 2 ), also depend on the gauge parameter. 
The gap equation's solution has the form 

Sfip)- 1 = i 1 -pA } (p\( 1 ) + B } {p 2 X 2 ). (5) 

and the mass function M f (p 2 ) = B f {p 2 ,C, 2 )/A f {p 2 X 2 ) 
is renormalisation point independent. The propagator is 
obtained from Eq. ([3|) augmented by a renormalisation 
condition. Since QCD is asymptotically free, the chiral 
limit is defined by 



Z 2 (( 2 ,A 2 )m^(A)=0, VA»C- 



(6) 



which is equivalent to requiring that the renormalisation 
point invariant current-quark mass is zero; i.e., rht = 0. 
A mass-independent renormalisation scheme can then be 
implemented by fixing all renormalisation constants in 
the chiral limit [18j |: namely, one requires 



SfipY 



lp2 = C 2 



l"f ■ P . 



(7) 



NB. 



Z 2 (C 2 ,A 2 )m^(A) = Z 4 (C 2 ,A 2 )m / (C). (8) 

While it is not readily apparent from Eqs. ©, (|3]), if 
one considers a Nf flavour theory, then there is coupling 
between the gap equations for different flavoured quarks. 
The coupling is driven by vertex corrections, as illus- 
trated with a model in Ref. Absent this coupling 
it is impossible to obtain other than a mean field chiral 
symmetry restoring phase transition at nonzero temper- 
ature [11]. We will return to this elsewhere but omit 
further substantial consideration herein. Consequently, 
we usually omit the flavour label in all that follows. 

The chiral susceptibility measures the response of a 
chiral order parameter to changes in current-quark mass. 
Consider therefore 



x(0 



d 



ct(to;C, A) 



dm(() 
Z 4 {(,A)N c tr D 



dm 



S(q;m;Q 



(9) 



(10) 



m=0 



The order of integration and differentiation can be inter- 
changed because the theory is properly regularised. One 
can now use the Ward identity 



d 



S(q; m; Q = ~S(q; to; ()T (q, 0; ()S(q; m; Q (11) 



dm 
to obtain 

X (0 = -Zt(C, A)iV c tr D f S(q; 0; ()r (q, 0; ()S(q; 0; C) , 

Jq 

wherein Tq is the renormalised fully-dressed scalar vertex, 
which satisfies an inhomogeneous Bethe-Salpeter equa- 
tion: 

r (fc,P;C) = Z 4 I D 



[S(q+)r (q, P)S(q-)} sr K™(q, k; P) . (13) 



In this equation: k is the relative and P the total momen- 
tum of the quark-antiquark pair; q± = q ± P/2; r, s, t, u 
represent colour and Dirac indices; and K is the fully- 
amputated quark-antiquark scattering matrix. 

In QCD the quantity m 2 ujo(P) is a renormalisation 
point invariant, where too(P) is the scalar vacuum polar- 
isation. In this product the term which multiplies m(C,) 2 
is 



1 In our Euclidean metric: {7 M ,7„} = 2<5 M „; 7^ = 7 M ; 75 
74717273; a ■ b = X)j=i a ibf, and P M timelike => P 2 < 0. 



uj (P;0 = Z 4 N c tv D S(q + )T (q,P)S(q 



(14) 
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Upon comparison with Eq. (| 1 2[) one arrives at the general 
result: 



x (C) = -a;o(P = 0;m = 0,C) 



(15) 



Hitherto we have not specified a regularisation proce- 
dure for the in-vacuum chiral susceptibility. In this con- 
nection it is noteworthy that if a hard cutoff is used, then 
in the chiral limit of a noninteracting theory 

-c °(F = 0;m = 0) = ^A 2 . (16) 

Following Ref. [17|, this result can be traced to the de- 
pendence on current-quark mass in Eq. On the 
other hand, Pauli-Villars regularisation would yield zero 
as the result, and this is the procedure we recom- 
mend and employ in models that preserve the one-loop 
rcnormalisation-group behaviour of QCD. 

One may implement a Pauli-Villars regularisation by 
introducing a pseudo-quark with large mass m p Q = A, 
which is anticoupled to the scalar source [35j]. To be 
concrete, this means 

w (P;to = 0,C) - lim Z 4 (C,A) 

A— »oo 



x N c tr D 



(2tt¥ 



S(q + ;0)T (q,P;0)S(q.;0) 



-S(q + ;A)T (q,P;A)S(q-;A) 



To proceed, one solves the gap equation, for the regular- 
quark and pseudo-quark, then uses the results to obtain 
the associated scalar vertices from Eq. ([13")) , and finally 
computes x{() from Eq. (fT5| . 

For the calculation of x(C) one onr y requires the scalar 
vertex at P = 0, at which total momentum it has the 
general form 



r (fc, 0; m) = ij ■ k C(k 2 ; m) + D(k 2 ; m) . 



(18) 



Owing to Eq. (jlip it is not necessary in principle to solve 
the inhomogeneous Bethe-Salpeter equation because 



C(k 2 ;m) 



_d_ 

dm 



A(k 2 ;m) ; D(k 2 ;m) 



_d_ 

dm 



B(k 2 ;m) 



Thus a solution of the gap equation suffices completely 
to fix Fo(k, 0). The practical utility of this procedure is 
limited only by the time required to solve the gap equa- 
tion numerically for a range of current-quark mass values 
and therefrom construct C, D. 

As another application of the scalar Ward identity, con- 
sider that in Landau gauge QCD one can write [2(J, HH 



A(k 2 ; m) 
B(k 2 -m 



1 . 



fe 2 »A| c 



(iln[fc7A3 QCD ])7„ 



{qqf 



fc2(|l n [fc2/A 2 CD ] 



(20) 



(21) 



where -f m — 12/ (33 — 2Nf), with Nf the number of active 
flavours, and (qq) is the renormalisation-group-invariant 
chiral-limit vacuum quark condensate. It therefore fol- 
lows from Eq. (Tl"Tj) that C = and 



D(k 2 ;m) k>> ^ QCD " ■ 2 <-> 



z 4 (C 2 ,fc 2 ) 



2^ 2 7r, 



where 



3 ^(|ln[feVA 2 CD ])i-7 m ' 
d(qq)° 



X Q = 

dm(Q 



(22) 



(23) 



Inserting Eqs. HID - (J23JI in Eq. (T7|) and using the fact 
that 37 m > 1 for all Nf, then 



-w o (P = 0;m = 0,C) = lim Z 4 ((,A)N C 



A 0-^2 



2^ 2 7« 



g2(| m [gVA 2 QCD ])i-^' 



(24) 



where all terms in the integrand that integrate to a finite 
value have been dropped. Thus 



(17) - MP = 0; m = 0, C) = A c ° ln[C/A QCD ]p 



d(qg) c 

dm (C) 
(25) 

This brief analysis demonstrates the consistency of our 
definitions and procedure. It also emphasises that the in- 
vacuum chiral susceptibility is a truly well-defined quan- 
tity. Indeed, one sees that in principle its value is con- 
tained in the solution for the dressed scalar vertex. 



III. GAP EQUATION MODELS 

A dialogue between DSE studies and results from nu- 
merical simulations of lattice-regularised QCD is provid- 
ing important information about the kernel of QCD's 
gap equation; e.g., Refs. [H, H, H H [H, E3, H HH . 
This body of work can be used to formulate reasonable 
Ansatze for the dressed-gluon propagator and dressed- 
quark-gluon vertex in Eq. (U) . 

In connection with such Ansatze it is notable that 
the DSEs admit at least one nonperturbative symmetry- 
preserving truncation scheme 0, [3(1 HH , which has en- 
abled the proof of numerous exact results [H, [H, [H, |H, 
HI, H3, HE HE EE EH . It also provides a starting point 
for the formulation of reliable models that can be used to 
illustrate those results and make predictions with readily 
quantifiable errors @, EE 51, EH EE HE S3 ■ The Ansatze 
are typically implemented by writing 



Z ig 2 D pa {p-q)T a a (q,p) 



G((p q) 2 ) D^ C (P q)^r a (q,p) , (26) 
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wherein D { ™ c (l) is the Landau-gauge free gauge-boson 
propagator, Q(l 2 ) is a model effective-interaction and 
T tT (q,p) is a vertex Ansatz. 

In one widely used approach to in- vacuum physics, 
G(£ 2 ) is chosen such that the one-loop renormalisation 
group behaviour of QCD is preserved and the vertex is 
written 



T<r(q,p) = la 



(27) 



This is the basis for a renormalisation-group-improved 
rainbow-ladder (RL) truncation of QCD's DSEs [H, EE 

One can alternatively employ Ansatze for the vertex 
whose diagrammatic content is unknown. A class of such 
models that has hitherto seen much use can be charac- 
terised by [48[ 



iT a {k, I) = iS A (k 2 ,£ 2 ) Ja + (k + i) 



'-^■(k + £) A A (k 2 ,e 2 ) + A B (k 2 ,e 2 ) 



where 



X F (k 2 ,£ 2 ) = ~[F{k 2 ) 



F(£ 2 )} 



A F (k 2 J 2 ) 



> , F(k 2 ) /••(/•-') 



k 2 - P 



(28) 

(29) 
(30) 



with F — A,B; viz., the scalar functions in Eq. ([5]). This 
Ansatz satisfies the vector Ward-Takahashi identity and 
is often referred to as the BC vertex. 

One has a Slavnov- Taylor identity for the quark-gluon 
vertex in QCD, not a Ward-Takahashi identity. Hence, 
Eq. (|28[) is not necessarily an improvement over Eq. ([2T| . 
A comparison between results obtained with the different 
Ansatze is nevertheless useful in identifying those out- 
comes which might be robust. 

Herein we employ a simplified form of the 
renormalisation-group-improved effective interaction 
in Refs. [HE!, H H EE E3; viz., we retain only that 
piece which expresses the long-range behaviour (s = k 2 ): 



Q{s) 



Att 2 



(31) 



This is a finite width representation of the form intro- 
duced in Ref. [4t|, which has been rendered as an inte- 
grable regularisation of l//c 4 [50] . Equation (j3"Tj) delivers 
an ultraviolet finite model gap equation. Hence, the reg- 
ularisation mass-scale can be removed to infinity and the 
renormalisation constants set equal to one. 

The active parameters in Eq. (|3ip are D and ui but they 
are not independent. In reconsidering a renormalisation- 
group-improved rainbow-ladder fit to a selection of 
ground state observables [Hj|, Ref. [IB)] noted that a 
change in D is compensated by an alteration of u). This 
feature has further been elucidated and exploited in 



TABLE I: Results obtained for selected quantities with u> = 
0.5 GeV, and the vertex and D parameter value indicated: 
A(0), M(0) are p = in- vacuum values of the scalar functions 
defined in connection with Eq. (JS| ; the vacuum quark conden- 
sate is defined with m = in Eq. ([2]); and \ is obtained from 
Eq. (|42|l . A(0) is dimensionless but all other entries are quoted 
in GeV. The calculations reported herein were performed in 
the chiral limit. 



Vertex 




A(0) 


M(0) 


-«W)°) 1/3 


Vx 


Eq. (27). RL 


i 


1.3 


0.40 


0.25 


0.39 


Eq.pg). BC 


i 

72 


1.1 


0.28 


0.26 


0.28 



Refs. [la, \M, |47J. For u G [0.3,0.5] GeV, with the in- 
teraction specified by Eqs. fl2BJ), J37J and (JHU), fitted in- 
vacuum low-energy observables are approximately con- 
stant along the trajectory 



luD = (0.8 GeV) s 



(32) 



Herein, we employ u> = 0.5 GeV. Therefore D — m z g /u> = 
1.0 GeV 2 corresponds to what might be called the real- 
world reference value for Eq. (|2T|) . 

It is impossible at present to explore whether the be- 
haviour expressed in connection with Eq. (|32p is also re- 
alised with the interaction specified by Eqs. (f2"o| . (28p 
and (PT|) because the BC vertex, Eq. (f2"8")) , cannot yet 
be used in the implementation of a symmetry preserving 
DSE truncation. We note, however, that D = 0.5 GeV 2 
reproduces the physical value of the in-vacuum conden- 
sate, last row of Table |U We therefore identify this as 
the real- world reference value for Eq. (|2"5|) . 

In association with Eqs. (fT9]) we noted that it might 
be time consuming to solve the gap equation numerically 
for a range of values of current-quark mass and therefrom 
construct C, D. It is therefore useful if that procedure 
can be avoided. In a class of cases, it can. 

From the gap equation 



r (p, o) = J^s-\ P ) = i D + ^-s(p) . 

om am 



One can always write 
d 



— E(p)=Lj(p,0) + L 2 (p,0), 



(33) 



(34) 



where, using Eq. pip . 

r5(p,0) = - / g 2 D fW (p-q) 

Jq 

rg(p,0)= J fD^ip-q) 

x y7^(«)y^r„( g , P ). 



(35) 



(36) 
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Now, with any Ansatz for the dressed-quark-gluon ver- 
tex whose current-quark-mass-dependence is completely 
determined by that of the dressed-quark propagator, 
■^T u (q,p) has a closed form involving C, D. In these 
cases Eq. (l33l) is a self-contained pair of coupled integral 
equations for C, D. NB. Eqs. p7|l. ([28]) are both in this 
class. In the former case Tq(p, 0) = 0, while in the latter 



_d_ 

dm 



-7 • (q + p) A c {q 2 ,p 2 ) + A D (q 2 ,p 2 ) 



(37) 



As noted above, the interaction in Eq. (|3Tj) does not 
preserve the one-loop behaviour of QCD but, indeed, 
delivers an ultraviolet finite gap equation. Neverthe- 
less, Eq. (Til?]) demonstrates that a regularisation proce- 
dure would be required for the vacuum polarisation even 
if one were dealing with a free field theory. For QCD, in 
which Schwingcr functions carry anomalous dimensions, 
the Pauli-Villars scheme works well. However, that pro- 
cedure, as defined above, fails when the vacuum polar- 
isation is evaluated with Schwinger functions generated 
by Eq. (|3"Tj) . A simple alternative can be constructed. 

The vector vacuum polarisation is given by 

,A 

^0P)=^ctr D / i-y ll S(q + )iT v (q,P)S(q-), (38) 

Jq 

where T„ is the vector quark-antiquark vertex; i.e., the 
vector analogue of r , which satisfies an inhomogeneous 
Bethe-Salpeter equation similar to Eq. (fT3|) . 2 This polar- 
isation couples to the photon and hence must be trans- 
verse and contain no mass term. Consider therefore [ill ] 

i at r A 

-co^(P = 0) = -^tr D J i^S^iT^q, 0)S(q) 

N °* f A - d q( ^ 

= -^jf^^Ms 8 )]. (39) 

where the vector Ward identity was used in the second 
line and we have written the dressed quark propagator 
in the form 



S(p) = -ij- P (T V (p 2 ) +a s (p 2 ). 



(40) 



Evaluating Eq. (|39[) with a hard cutoff in a noninteracting 
theory one obtains 



(41) 



i.e., precisely the same sort of divergence as encountered 
in the scalar vacuum polarisation. Naturally, any regular- 
isation scheme that preserves the vector Ward-Takahashi 
would yield zero as the result [5l| . 

In connection with Eq. (f3"Tj) we therefore define the reg- 
ularised in- vacuum chiral susceptibility as 

X = -wo(0; m = 0, A) + ^(0; m = 0, A) (42) 

at r°° ( 
= --^j dsslD(s)[a s (s) 2 -sa s (s) 2 } 

+2s C(s)a v (s)a s (s) + 2a v (s) + sa' v (s) \ (43) 



2 Recall that with Eq. 1311 the regularisation mass-scale should be 
removed to infinity and the renormalisation constants set equal 
to one. 



It will be recognised that so long as each term in Eq. (j42|) 
is independently regularised in a valid fashion, which one 
can safely assume, then this is a nugatory transforma- 
tion. After all, we have only subtracted zero: the pho- 
ton is massless and hence in any valid scheme it must 
be that (0;0,A) = 0. On the other hand, with this 
definition one can proceed to calculate the susceptibility 
without it being necessary to actually specify a regular- 
isation scheme. (NB. The last two contributions in the 
integrand act as the Pauli-Villars terms. Z4 = 1 in the 
ultraviolet-finite model.) 



IV. ILLUSTRATION 

It is now a straightforward matter to solve for the 
dressed-quark propagator and the P = scalar vertex. 
Owing to the GauBian form of Eq. (f3"Tj) , all relevant inte- 
grations converge rapidly. In Fig. [1] we plot the functions 
obtained through solving the gap equation and in Fig. [2] 
those which describe the P = scalar vertex. The results 
were obtained with the appropriate real-world value of 
the interaction strength: D = 1 GeV 2 with Eq. (f2"7|) and 
D = 0.5 GeV 2 with Eq. 

It is apparent in Fig.[T] that the vertex Ansatz has a 
quantitative impact on the magnitude and pointwise evo- 
lution of the gap equation's solution. That this should 
be anticipated is plain from Ref . [53] ■ Moreover, the pat- 
tern of behaviour can be understood from Ref. [24j : the 
feedback arising through the A# term in the BC ver- 
tex, Eq. |25]) , absent in Eq. JS7J), always acts to alter the 
domain upon which A(p 2 ) and M(p 2 ) differ significantly 
in magnitude from their respective free-particle values. 
Since C(p 2 ) and D(p 2 ) are derived quantities, their be- 
haviour does not require explanation. 

The integrand in Eq. (|4"3"]l is depicted in Fig. [3] for each 
vertex Ansatz at the associated real-world interaction 
strength. The resulting chiral susceptibilities are pre- 
sented in Table[l] Two things are immediately apparent 
in the figure. First, the regularisation has served to elim- 
inate the far-ultraviolet tail of the integrand, thereby en- 
suring convergence of the integral. We have varied the 
detailed form of the regularising subtraction; e.g., us- 
ing free-field propagators and vertices instead of gap and 
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FIG. 1: Dressed quark propagator. Upper panel - A(p ) = 
l/Z(p 2 ), where Z is the wave function renormalisation func- 
tion. Lower panel - B(p 2 ), the scalar piece of the dressed- 
quark self-energy. In both panels, Dashed curve : calculated 
in rainbow-ladder truncation, Eq. (|27)l. with I — D/lj 2 — 4; 
solid curve : calculated with BC vertex Ansatz, Eq. (|28[) , and 
X = 2. 
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Q 




10"" 10" 3 10" 2 10" 1 10° 10 1 10 2 

p 2 (GeV 2 ) 

FIG. 2: P — scalar vertex, Eq. (flHl upper panel - C(p 2 ), 
lower panel - D(p 2 ). In both panels, Dashed curve : cal- 
culated in rainbow-ladder truncation, Eq. (|27[) , with I = 
D/lj 2 = 4; solid curve : calculated with BC vertex Ansatz, 
Eq. |(SH), and 1 = 2. 



Bethe-Salpeter equation solutions. The pointwise be- 
haviour of the integrand is little altered. Second, at real- 
world values of the interaction strength for both Ansatze 
the integrand has negative support in the infrared and 
positive support for p 2 > (0.5 GeV) 2 . We will subse- 
quently return to this point. 

In Fig. 2] we depict the evolution of the chiral suscep- 
tibility with increasing interaction strength, X — D/lo 2 . 
The behaviour may readily be understood. For 1 = one 
has a noninteracting theory and the "vacuum" is unper- 
turbed by a small change in current-quark mass. Hence, 
the susceptibility is zero. The susceptibility grows with 
increasing X because the interaction is attractive in the 
qq channel and therefore magnifies the associated pair- 
ing. This is equivalent to stating that the scalar vertex 
is enhanced above its free field value. 

The growth continues and accelerates until, at some 
critical value, X = X c , the susceptibility becomes infinite. 



Those critical values are: 

Eq. (EH), RL 
l r 1.93 



Eq. dMJ), BC 
1.41 



(44) 



This divergence is easily understood. The models we've 
defined contain a dimcnsionless parameter, X, which 
characterises the interaction strength, and a current- 
quark mass, which is an explicit source of chiral symme- 
try breaking. In the general theory of phase transitions 
the latter is analogous to an external magnetic field while 
1 jX is kindred to a temperature. 

Consider the free energy for such theories f(t,m), 
where t = [1—X C /X]. If such a theory possesses a second- 
order phase transition, then the free energy is an homo- 
geneous function of its arguments in the neighbourhood 
of t = = m. From this it follows that the theory's 
magnetisation exhibits the following behaviour (e.g., see 
the appendix of Ref . [l(| ) : 

M(t,0) =t ,t^ 0+ (45) 



p 2 (GeV 2 ) 

FIG. 3: Integrand in Eq. (|43[) - Dashed curve : calculated in 
rainbow- ladder truncation, Eq. (127[) , with X = D/lj 2 = 4; 
solid curve : calculated with BC vertex Ansatz, Eq. (|28[) , and 
1 = 2. 



and the associated magnetic susceptibility evolves ac- 
cording to 

M(0,m) oc m- {1 - 1/s) . (46) 

NB. For a mean- field theory (3=1/2 and <5 = 3. 

The nature of the critical interaction strength is now 
plain. In the class of theories we're considering, the vac- 
uum quark condensate is analogous to the magnetisa- 
tion. It is attended by the chiral susceptibility. 3 For 
I <I C the interaction has insufficient strength to gener- 
ate a nonzero scalar term in the dressed-quark self-energy 
in the absence of a current-quark mass; namely, dynam- 
ical chiral symmetry breaking (DCSB) is impossible and 
the model realises chiral symmetry in the Wigner-Weyl 
mode. That changes at X c , so that for X > X c a B ^ so- 
lution is always possible. Moreover, the behaviour of the 
susceptibility shows that each model undergoes a second- 
order phase transition and realises chiral symmetry in the 
Nambu-Goldstonc mode for interaction strengths above 
their respective values of X c . 

These observations emphasise the chiral susceptibil- 
ity's usefulness. With the vertex of Eq. (f2"T)) one can ex- 
plicitly construct the pressure and therewith show that 
for X > X c the DCSB solution is dynamically favoured 
because it corresponds to the configuration of maximum 
pressure. On the other hand, the diagrammatic content 



FIG. 4: Dependence of the chiral susceptibility on the inter- 
action strength in Eq. (|31[l ; viz., X := D/lo 2 : dashed curve, 
RL vertex, Eq. (J2JJ; solid curve, BC vertex, Eq. (25) . 



of the vertex in Eq. (|28|) is unknowable and hence one can- 
not derive an expression for the pressure. In this case one 
may rely on the behaviour of the susceptibility to con- 
clude that DCSB is favoured, as illustrated in Ref. [54| 
within the Nambu-Jona-Lasinio model. 

Once past the critical point, the susceptibility de- 
creases as X increases. This is because the magnitude 
of the condensate order-parameter grows in tandem with 
X and therefore the influence of any perturbation associ- 
ated with a current-quark mass must steadily diminish. 5 

As in the earlier figures, Fig. 2] displays results with a 
quantitative dependence on the vertex Ansatz. This is 
easily understood [241 ]: at any given, common interac- 
tion strength in Eq. (|3"Tj) . the BC vertex Ansatz, Eq. (|28j) . 
amplifies attraction in the qq channel in comparison with 
that obtained in the rainbow truncation, Eq. (127[) . 

In Fig. © we depict the evolution with X of the inte- 
grand in Eq. (|42|) . This illustrates the manner by which 
the susceptibility's behaviour is generated. For X < X c 
the integrand is a nonnegative monotonically decreas- 
ing function whose pointwise magnitude grows rapidly 
as X — » Plainly, given the scalar Ward identity, 

Eq. (fTTj) . the functions characterising the scalar vertex 
are the origin of this rapid growth: the vertex exhibits 
a singularity at X c . For X > X c the integrand possesses 
a zero. Its origin is the existence on this domain of a 
nonzero mass function in the chiral limit, which provides 
for a negative contribution owing to the D cr| term in 

The behaviour of the scalar vertex and the susceptibil- 
ity that we have described contains additional, valuable 
information. At X = one has a free field theory and 



3 The dressed-quark mass function evaluated at p 2 = and the 
pion's leptonic decay constant are equivalent order parameters. 

4 The pressure is defined as the negative of the effective-action. A 
system's ground state is that configuration for which the pressure 
is a global maximum or, equivalently, the effective-action is a 
global minimum. These statements are elucidated in Ref. |53|| . 



5 On the other hand, if at fixed coupling one were to increase 
the current-quark mass from zero, the dynamical component of 
the order parameter would steadily be suppressed in response 
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FIG. 5: Influence of X := D/lo 2 on the pointwise behaviour 
of the dimensionless integrand in Eq. (|43[) . evaluated with the 
BC vertex, Eq. (128]! . 



no bound state in the scalar channel, nor any other. As 
X increases away from zero, attraction appears in the 
scalar channel and with it a qq correlation. Nevertheless, 
for X < X c in the chiral limit the dressed-quarks remain 
massless. Hence, absent confinement, zero is the only 
realisable mass for a bound state in the scalar channel. 
However, since \ involves the scalar vertex evaluated at 
P = and < x < 00 j it is apparent that no scalar 
bound state exists when X < X c . 

On the other hand, the chiral limit susceptibility ex- 
hibits a singularity at X = X c . This feature is tied to a 
pole at P 2 = that appears in the scalar vertex. Thus at 
the critical interaction strength a massless bound state 
appears in the scalar channel. Moreover, since at pre- 
cisely X c the dressed-quark mass is still zero, then a de- 
generate; i.e., massless, pseudoscalar bound state appears 
simultaneously It is evident in the evolution of \ that 
the scalar bound state becomes steadily more massive as 
X increases away from the critical value. 

In terms of t = [1—X C /X] the pattern of behaviour just 
described is similar in many respects to that exhibited 
by such models at nonzero temperature; e.g., Ref. [57l ]. 
The key difference is the absence of anything with the 
character of a screening mass for X < X c . The existence 
of a screening mass allows a pole in the scalar vertex for 
temperatures above critical. 



V. SUMMARY 

We demonstrated that the in-vacuum chiral suscepti- 
bility is given by the scalar vacuum polarisation at zero 
total momentum. The scalar Ward identity is essen- 
tial in deriving this result from the expression for the 
vacuum quark condensate. One can rigorously define 
the in-vacuum susceptibility by employing a Pauli-Villars 



regularisation procedure. While this provides the most 
practical means by which to calculate the susceptibility, 
the derivation shows that, in principle, all information 
about the susceptibility is contained within the P = 
behaviour of the dressed scalar vertex. This confirms a 
posteriori that a valid regularisation of the scalar vacuum 
polarisation is always possible. 

One must solve the gap equation in order to calculate 
the chiral susceptibility. The kernel of that equation in- 
volves the dressed quark-gluon vertex. If an Ansatz for 
that vertex is employed whose diagrammatic content is 
unknown, then a consistent Bethe-Salpeter kernel cannot 
currently be constructed. In general this prevents the 
calculation of all vertices and in particular those that 
couple to colour singlet q-q channels. However, owing 
to the Ward identities, that problem is circumvented at 
P = in the scalar and vector vertices. This is a mate- 
rial point. For example, herein it means that no mat- 
ter what form is assumed for the dressed-quark-gluon 
vertex and, indeed, for the kernel of the gap equation, 
the consistent chiral susceptibility can always be con- 
structed. Additional simplifications follow if one employs 
an Ansatz for the quark-gluon vertex whose current- 
quark-mass-dependence is completely determined by that 
of the dressed-quark propagator. These features enabled 
us to provide a consistent calculation of the chiral sus- 
ceptibility for two rather different vertex Ansatze. The 
results were consistent and the minor quantitative differ- 
ences easily understood. 

A key result is that the susceptibility can be used 
to demarcate the domain of coupling strength within a 
theory upon which chiral symmetry is dynamically bro- 
ken. For couplings below the associated critical value 
and in the absence of confinement, there are no bound 
states. Degenerate massless scalar and pseudoscalar 
bound-states appear at the critical coupling. The be- 
haviour of the susceptibility shows that the scalar bound- 
state becomes massive as the coupling increases above 
this value whereas, naturally, the Goldstone mode re- 
mains massless in the chiral limit. 

We described and used a simple model for the gluon 
propagator in the gap equation and hence in computing 
the susceptibility. Nevertheless, no qualitative feature 
of our results will change should a form be used that, 
for example, expresses properly the one-loop renormal- 
isation group behaviour of QCD. That is demonstrated 
by calculations which are now underway. 

In addition, effects that might arise through flavour 
mixing in the gap equation are quashed by the vertex 
Ansatze we considered. Such mixing is generated by non- 
trivial correlations in the vertex. It is certainly impor- 
tant in-medium and should also be considered in-vacuum, 
where it is likely to alter the behaviour of the susceptibil- 
ity in the neighbourhood of the critical coupling strength. 
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